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Abstract 

Double-7 vibrations in deformed nuclei are analyzed in the context of the 
interacting boson model. A simple extension of the original version of the 
model towards higher-order interactions is required to explain the observed 
anharmonicities of nuclear vibrations. The influence of three- and four-body 
interactions on the moments of inertia of ground- and 7-bands, and on the 
relative position of single-7 and double-7 bands is studied in detail. As an 
example of a realistic calculation, spectra and transitions of the highly 7- 
anharmonic nuclei ^^^Dy, ^^^Er, and ^^^Er are interpreted in this approach. 
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I. INTRODUCTION 



Vibrational degrees of freedom in atomic nuclei can be described in terms of phonon 
excitations that arise from nuclear shape oscillations. Vibrations of nuclei with ellipsoidal 
symmetry can be of two types |l|| : (3 vibrations which preserve axial symmetry and give rise 
to a band with K = 0, and 7 vibrations which break axial symmetry and yield a. K = 2 
band, where K is the projection of the angular momentum on the axis of symmetry. At the 
experimental level, 7 bands have been identified in many well-deformed nuclei; in contrast, 
the identification of (3 bands is still full of questions and difficulties. This is mainly because, 
when the energy surface has a well-deformed minimum in /3 but is rather flat in 7, the /3- 
band increases in excitation energy and approaches the energy region where other degrees of 
freedom are important. In that case band mixing may occur and can give rise to non-pure 
structures with decay patterns difficult to identify as those of a /5 band 0. 

Since single-7 excitations are very well established, it is natural to look for double-7 vi- 
brations and to develop models that can deal with such multi-phonon excitations. Double-7 
excitations correspond to K'^ = 0^ and = 4+ bands which are the anti-parallel and paral- 
lel combinations of single-7 phonons, respectively. The experimental identification of two-7 
states in deformed nuclei is difficult because their expected excitation energy is around the 
pairing gap and hence they can mix strongly with two-quasiparticle excitations. However, 
recent experimental improvements in nuclear spectroscopy following Coulomb excitation , 
inelastic neutron scattering [Q, and thermal-neutron capture have made possible the 
study of highly excited low-spin states. Many states have been proposed as possible candi- 
dates of double-7 vibrations. There is, however, some controversy about their interpretation. 
The author of reference claims that some of the presumed double-7 states can be inter- 
preted as single hexadecapole-phonon excitations. In fact, to identify the band-head of a 
double-7 band it is not sufficient to analyze just B{E2) values; data from single- nucleon 
transfer reactions, /?-decay studies, and inelastic scattering experiments must be considered 
as well. One of the key properties to disregard a band as a double-7 band is the fact that 
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its members, in first order, cannot be populated in single-nucleon transfer reactions. Many 
examples of K'^ = states tliat are identified as double-7 excitations but are strongly 
populated in single-nucleon transfer reactions, can be found in the literature: ^^'^Gd, ^^^Dy, 
^^^Yb, ^''^'^^'^Hf, and ^^'^'^^^Os. Since the double-phonon character of the states in question 
is in doubt, they are not considered here. However, some candidates seem to have a gen- 
uine double-phonon nature. Such is the case with ^^^Dy and i^e-esjr^j, particular, in 
reference a = 4"'" state in ^^^Dy at 2.173 MeV is found to exhibit all properties of 
a double-7 band. In references P,^] the observation is reported of K'^ = O"*" and = 4"'" 
double-7 states in ^^^Er, at energies of 1.949 MeV and 2.029 MeV, respectively. Finally, in 
^^^Er a K'^ = 4+ double-7 excitation is identified at an energy of 2.055 MeV [p!0[] . 

One of the most striking features of the observed double-7 bands is their high anhar- 
monicity, i.e. the ratio of double-7 over single-7 energy is different from 2 and ranges from 
2.5 to 2.8. This information is very important since it provides a stringent test of nuclear 
models. The nuclei ^^"^Dy and i^e.iesgj, have been interpreted in the context of many different 
models such as the quasi-phonon model [ PT| , P^ , the geometrical model |T^, the multi-phonon 
model |T^, the self-consistent collective-coordinate method |T^|T^, and the sdg-lBM []T^,|TB[, 
and it is now of interest to revisit these models in connection with anharmonic vibrational 
behavior. 

In this paper the simplest version of the interacting boson model (IBM) |jT^ is extended 
by adding to the usual Hamiltonian higher-order interactions between the bosons with the 
purpose of creating a framework that accommodates the high anharmonicities observed in 
^^^Dy and i66,i68j^j,_ rjj^g structure of the paper is as follows. First, the IBM is reviewed 



with special reference to its harmonic character. In section |ITl| the inclusion of three-body 
terms in the Hamiltonian is discussed. The introduction of four-body terms is presented 
in section |rV| and some analytic results are pointed out. In section |V| a detailed study of 
possible four-body terms is carried out and realistic calculations for ^^^Dy and i66,i68gj, 
presented. Finally, in section the conclusions of this work are made. 
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II. THE IBM-1 AS A HARMONIC MODEL 



The IBM describes low-lying collective excitations in even-even nuclei in terms of 
monopole (s) and quadrupole (d) bosons [IS]. The boson number that corresponds to a 
given nucleus equals half the number of valence nucleons (A^ = n/2). The rotationally 
invariant and number- conserving boson Hamiltonian usually includes up to two-body in- 
teractions between the bosons although higher-order terms can be added in principle. The 
most general two-body IBM Hamiltonian can be written in a multipole expansion as 

H = Eshs + Sdhd + kqP^ P + KiL ■ L + K2Q ■ Q + K3T3 ■ T3 + K^f^ ■ T4, (1) 

where fig and fid are the s- and d-boson number operators, respectively, and 

Pt = irft . rft _ l,t . (2) 
L = yT0(rft X (3) 
Q = stj+ d^S + x{d^ X J) (2), (4) 

Tg = (rft X df\ (5) 

f4 = (rft X (6) 

The symbol ■ represents the scalar product; in this paper the scalar product of two oper- 
ators with angular momentum L is defined as Tl ■ Tl = J^Aii''^)^^ T^mTl-m where Tlm 
corresponds to the M component of the operator Tl- In the previous equations the operator 
Jem = (— l)™7f-m (where 7 refers to s or d) is introduced so that the annihilation operator 
verifies the appropriate properties under spatial rotations. 

It is not a priori clear to what extend f3 and 7 vibrations are anharmonic in the IBM 
even if one just considers the Hamiltonian (|T]) with up to two-body interactions. A partial 
analysis of this problem was given in reference |^ . There, the authors find that the IBM in 
its simplest version is a harmonic model in the limit of infinite boson number and even for 
finite the model cannot accommodate large anharmonicity if one considers up to two-body 
interactions; only the interplay between one-|-two-body terms and higher-order interactions 



can induce, in principle, a sizeable anharmonicity in the double-phonon excitations. The 
reason why one-body terms and two-body interactions cannot create a large anharmonicity 
can be understood as follows. If one considers a Hamiltonian with one parameter that 
controls the ratio of the strength of the one-body energies and the two-body interactions, 

where ^ ranges from to 1, one finds two 'phases' separated by a critical value, a first 
phase where the one-body term plays the main role < ^c) and a second phase where the 
two-body interaction is the driving force > .^c)- The crucial point is that the separation 
between the two phases is very sharp and essentially no interplay between one- and 



two-body terms can be found. Since, to a good approximation, the force is either one body 
or two body but not both, harmonic behavior cannot be avoided. 

The inclusion of high-order interactions in a system with a high boson number also leads 
to a harmonic description. Only for finite boson number the interplay between one-|-two- 
body terms and higher-order interactions can induce an anharmonicity in the double-phonon 
excitations that is comparable to the observed one. These ideas will be used as a guideline 
in the following sections. 

To carry out a quantitative study of anharmonicities, it is convenient to define a ratio 
between single- and double-phonon excitation energies. Because the experimental situation 
for (5 excitations is not clear we concentrate on 7 vibrations and define energy ratios for 7 
phonons only (although similar definitions can be given for j3 phonons): 



where 0+^ and 4+^ are the band heads of the K'^ = 0"*" and K'^ = 4+ double-7 bands, 
respectively, and stands for excitation energy. This particular definition removes any 
rotational influence. 
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III. THREE-BODY HAMILTONIANS 



Let us consider in the following a Hamiltonian that includes a quadrupole term, a rota- 
tional term, and three-body interactions between the d bosons, 

H = -kQ ■Q + k'L- L + Y^di ({d) X X d^Y^ ■ ({d X df^^ X d^^ , (9) 

kl 

where — k = K2, = i^i- Five independent three-body ci-boson interactions exist which 
have / = 0, 2, 3, 4, and 6. Interactions with the same I but different k are not independent 
but differ by a normalization factor only [^. The combinations (A;,/) = (2,0), (0,2), (2,3), 
(2,4), and (4,6) are chosen here. 

The Hamiltonian (^ is certainly not the most general one-|-two+three-body Hamiltonian 
that can be considered. Notably, a vibrational term edfid which dominates in spherical nuclei 
is omitted since it is thought of lesser importance in the deformed nuclei considered here. 
And, of all possible three-body interactions (seventeen terms), only those between the d 
bosons are retained here since they are the more efficient terms to produce anharmonicity 

For the discussion of the anharmonicities of 7 vibrations we study the behavior of the 
energy ratios (||) as a function of the ratio 61 j k. The identification of the states O;];^ and 4+^ 
is based on the B{E2) values for decay into the single-7 states. In figure ^ the infiuence 
of the various three-body interactions is shown for a typical value of x (x = —0.5) and for 
= 15 bosons. It is seen that a 7-vibrational anharmonic behavior is obtained which can 
be different for the = 0^ and = 4^^ bands {e.g. positive for the former while negative 
for the latter.) Care has been taken to plot results only up to values of 61 that do not 
drastically alter the character of rotational spectrum; beyond these values, the three-body 
interaction, being of highest order in the Hamiltonian (^, becomes dominant. Also shown 



^Note that this figure differs from figure 2 in reference |24| in some scale factors due to an error 
in the definition of 9i. 



in figure |l| are the ratios R]^ as observed in ^^^Er Rq = 2.76 and Rj = 2.50. Figure |^ 
shows the experimental spectrum of ^^^Er and compares it to the eigenspectrum of 
Hamiltonian (|^) with an Z = 4 three-body interaction. The parameters are n = 23.8 keV, 
X = —0.55, k' = —1.9 keV, and 64 = 31.3 keVQ, with boson number N = 15. With these 
values the calculated excitation energies of the double-7 band heads are 1926 keV and 1972 
keV for the K'" = 0^ and K'^ = 4"*" levels, respectively, leading to the ratios Rq = 2.82 
and Rj = 2.45, in excellent agreement with observation. Note, however, that although all 
7-band heads are well reproduced by the calculation, problems arise for the moments of 
inertia, in particular of the 7 band. In next sections we will come back on the moments of 
inertia to see that three-body Hamiltonians provide a very poor description of them. 

IV. SU{3) HAMILTONIANS WITH UP TO FOUR-BODY INTERACTIONS 



In the previous section and in reference |24| a very good description of double-phonon 



excitation energy has been obtained, but at expense of spoiling the moments of inertia of 
ground and 7 bands. These drawbacks seem to be a general feature of three-body Hamil- 
tonians. In this and the following section it is shown that the drawbacks of three-body 
interactions can be overcome by going to the next order. 

Since 7 anharmonicity has been observed exclusively in well-deformed nuclei, it is appro- 
priate to consider the problem in the SU(3) limit of the IBM which is suited to deal with 
nuclei in this mass region p5[. Therefore, in a first approach, the SU{3) limit is used and 
later, in the next section, these results are used as a guidance in more realistic calculations. 

Let us consider the following Hamiltonian: 

H = a C2[SUi3)] + 61 Cs[SUi3)] + 62 NC2[SUi3)] 

+ ci C2[SUi3)f + C2 NC3[SUi3)] + C3 N^C2[SUi3)], (10) 



^ Note that the value 64 is 1/3 of that given in reference |24] which has an error in its definition. 
The results shown in that paper are correct after a simple rescaling of the parameter. 
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where Cn[SU{3)] stands for the Casimir operator of order n of SU{3). Note the inclusion 
of cubic terms for completeness in the SU (3) analysis. Since only the spectrum of a single 
nucleus is of interest here, the number of bosons can be fixed in every case and the Hamil- 
tonian ([T0|) can be simplified by combining terms into a single one, leaving a Hamiltonian 
with three coefficients a, b, and c, 

H = a C2[SU{3)] + b Cs[SU{3)] + c C2[SUi3)]\ (11) 

with eigenvalues 

((A,/i)|if|(A,/i)) = a(A2 + + A/i + 3A + 3/i) 

+ 6(A - /i)(2A + ^ + 3)(A + 2/i + 3) 

+ c{X^ + fi^ + Xfi + 3X + 3fi)\ (12) 

No quartic Casimir operator exists for SU{3) because the number of independent Casimir 
operators equals the number of labels that characterize an irreducible representation. The 
Hamiltonian ([Tl|) has no rotational term since of primary interest, at this point, is the 
description of band-head energies of single- and double-7 excitations. 

The definition of the energy ratios (^ must now be adapted to incorporate the symmetry 
labeling of the states. In the SU (3) limit the 7 band belongs to the SU (3) representation 
{2N — 4, 2) where N is the number of bosons. The double-7 band with = 4 is contained 
in the (2A^ — 8,4) representation; the double-7 band with K = is predominantly contained 
in the (2A^ — 6, 0) representation, although an important component is in (2A^ — 8, 4) pO| . 
In this section the energy ratios (P) are thus defined as follows: 

5c/{3)7 ^ ^x(2iV-6,0) sr/(3)7 ^ ^x(2iV-8,4) 

° " E,(2iV-4,2)' ^ " E,(2iV-4,2)- ^ ^ 



Because no rotational term is included, the ratios (|T3D can be compared directly with equa- 
tions (^. In the following the effect of the different terms in equation ([TT|) on the degree of 
anharmonicity of the two-phonon excitation energy is analyzed. 
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a < 0, 6 = 0, c = 0. 



This corresponds to the simplest version of IBM; the energy ratios become 

sui3), 24iV - 18 su(3), 2AN - 36 

^° - 12iV-6' - 12N-Q- ^^^^ 

An almost pure harmonic 7-vibrational spectrum is found since the energy ratios (|n|) 
are only slightly lower than 2. This will be referred to as negative anharmonicity as 
opposed to the positive anharmonicity for energy ratios above 2. 

a<0,b^O,c = 0. 

In this case the Hamiltonian (|llj) is a combination of two- and three-body terms. For 
given values of a and b and for a high enough boson number, only the three-body part of 
the Hamiltonian is dominant and a harmonic spectrum is recovered. For obtaining an 
anharmonic spectrum the values of the Hamiltonian parameters are very constrained 
once the number of bosons has been fixed, as the following analysis shows. 

Without lost of generality a can be fixed to a = — 1. The energy ratios are then 

sui3), ^ 3-4iV + fe(24iV^-36iV + 27) 
° {2N-l){6bN + 9b-l) ' 

R. - 2^^^. (16) 

To keep the energy of the 7 excitation positive, the value of b has an upper limit 

b < . (17) 

The value b = 1/{6N + 9) leads to a divergence in Rq^^^''"' and around this value 
anharmonic behavior is found. From equation ([T6|) one observes that the behavior of 
the (2A^ — 8, 4) representation is completely harmonic and does not depend on b. On 
the other hand, from equation (|l^) one sees that a wide range of anharmonic ratios is 
found for the (2A^ — 6, 0) representation. As an illustration, in figure |^ equation (|15D 
is represented as a function of b, for three values of (5, 10, and 15). Only positive 



values of b are plotted because for the negative ones Rq'^^^^"' decreases smoothly to the 
asymptotic values 1.27, 1.58, and 1.70 for N = 5, 10, and 15, respectively. 

The conclusion is that a Hamiltonian with c = does not agree with the experimental 
situation observed in the mass region of well-deformed nuclei, where the 7|-=4 state is 
highly anharmonic. 

a < 0, 6 = 0, 0. 

In this case the Hamiltonian ([ll|) is a combination of two- and four-body terms. As 
in the previous case, anharmonicity requires very constrained Hamiltonian parameters 
once the number of bosons is fixed. 

Again, without loss of generality, we fix a = —1. The energy ratios then read as 
follows: 

5c/(3h ^ (4iV-3)(2c(4iV^-6Ar + 9)-l) 

° (2iV-l)(8cAr2 + 6c-l) ' ^ ^ 

p5c/(3h ^ 2(2iV-3)(4c(2Ar2-3iV + 9)-l) 
' {2N - l){8cN^ + 6c - 1) ■ ^ ^ 

To keep the energy of the 7 excitation positive, the value of c has the upper limit 

c < 1 . (20) 

The value c = 1/(8A^^ -|- 6) produces a divergence in the two energy ratios and in its 
neighborhood highly anharmonic behavior is found. In figures ^ and ^ are plotted the 
ratios R^'^^^^"' and R^'^^^^"' , respectively. Again, for negatives values of c, i?Q^''^^^ goes 
asymptotically to the values 1.45, 1.69, and 1.78 for = 5, 10, and 15, respectively, 
while goes to 1.33, 1.59, and 1.71. These negative anharmonicities have no 

phenomenological interest. 

In this case the experimental situation can be nicely described. Both Jx^q and 7|:=4 
states can be accommodated in a anharmonic description. The conclusion is thus that 
a Hamiltonian with C[SU{3)] and C'[S'[/(3)]^ terms seems to be a good starting point 
to treat the 7 anharmonicity in deformed nuclei. 
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The description presented here only provides band heads and, to recover a rotational 
structure, a U' term must be included. The rotational structure is the same in every band 
because in the SU{?)) limit no mixing exists between rotational and vibrational degrees of 
freedom. In next subsection these results are illustrated with a schematic calculation for 
energies and transition probabilities. 



A. A schematic application 

Let us consider the case of ^^^Er which is, as already mentioned, one of particular interest 
because both double-7 excitations (with K = Q and K = A) have been identified 
To carry out the schematic calculation, we use the Hamiltonian (|ll|) with 6 = 0. The 
experimental values for the single and double-7 energy ratios are i?Q = 2.76 and i?] = 2.50 
and can be compared directly with the expressions (|18|-|19|) leading two values for — c/a. 



namely 5 x 10~^ and 1 x 10^^. Both solutions are fairly close and any value in between them 
will correctly describe the anharmonicity of the K = {] and K = A bands. The value of a 
and the strength of the rotational term are fixed from the excitation energies of the 2i and 
22 levels. After these simple considerations one arrives at the following Hamiltonian: 

H = 13.43 L-L- 20.84 C2[SU{3)] + 9.296 10"^ C'2[5f/(3)]2, (21) 

where the coefficients are given in keV. The theoretical and experimental spectra are com- 
pared in figure ^ and a very good agreement is obtained. However, due to the simplicity of 
the calculation, 7 and j3 bands are degenerate in energy which is not the case experimentally. 

To complete the description, E2 transition probabilities must be computed also. The 
calculation of B{E2) values in the SU{3) limit, as in other situations where degeneracies 
occur, must be treated with care and an appropriate basis must be chosen for states with 
the same energy. A natural way to do this work is to slightly lift the degeneracy of the 
SU{3) Hamiltonian. The levels can be split using in the Hamiltonian a value x = —1.30 
which is very close to its SU (3) value. With this SU (3) breaking the degeneracy is lifted in a 
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natural way because the (3 band is pushed up in energy, as is observed. One may expect that 
with this small change in x the SU (3) spectrum will keep its properties. The E2 transition 
operator is: 



operator and in the Hamiltonian is different. The effective charge is fixed to reproduce 
B{E2;2f of): = (1.97)^ W.u.. In table | theoretical and experimental transition 
rates involving the ground and the 7 bands are compared. 

This simple analysis suggests that a four-body operator of the type in (^) provides a 
good description of both single- and double-7 bands. In the next section this schematic 
analysis is extended to non-SU (3) situations. 

V. GENERAL HAMILTONIANS WITH UP TO FOUR-BODY INTERACTIONS 

A general Hamiltonian with all possible three- and four-body terms can, in principle, 
be constructed but the number of parameters is so high that a study, even a schematic 
one, of the effect of the different terms on energy spectra and electromagnetic transitions is 
impossible. Schematic IBM Hamiltonians have been used for many years and, in particular, 
the quadrupole-quadrupole interaction has been very successful in describing a wide variety 
of nuclear spectra [1T9| , ^ . On the other hand, in the previous section it was shown that an 
expansion in terms of Casimir operators, which are mainly related to quadrupole operators, 
leads to a satisfactory description of ground, single- and double-7 bands. It is worth noting 
that such an expansion in terms of Casimir operators has been successfully used in Molecular 
Physics where spectroscopic data provide many anharmonic states ||2^. This leads us to 
propose a Hamiltonian as a quadrupole expansion that includes up to four-body terms. An 
alternative Hamiltonian can be based on an expansion in terms of pseudo-Casimir operators, 
which we define here as operators that become a true Casimir operator only for a particular 




(22) 



The value of x that best reproduces the data is x 



0.26. Note that x in the T{E2) 
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choice of one structure parameter. For example, the C*2[S'f/(3)] operator is related to the 
quadrupole operator through 2Q ■ Q + where Q = s'^d + Ss ± y/7/2{d^ x d)^'^\ If the 
value ±V^/2 is changed to a new operator C2[SU{3)1^ is obtained which we refer to as 
a pseudo-Casimir operator. It is a Casimir operator only for x = /2. 

Guided by the results of the previous section, two possible Hamiltonians that include up 
to four-body interactions, can be proposed, one based on a quadrupole expansion and the 
other on a pseudo-Casimir expansion: 

HQ = K'L-L + aQ-Q + biQxQx g)^ + c {Q ■ Q){Q ■ Q), (23) 
HpC = K' L-L + a C2[SU{3)]^ + b C^[SU{3)]^ + c C2[SU{3)]1, (24) 

where Q = Q{x) ^.nd 

C2[SU{3)]^ = 2Q{x)-Q{x) + Il\ (25) 

C3[SU{3)]^ = -4V35{Qix) x Qix) x Qixif'^ - ^v^(L x L x Qix)Y'\ (26) 

For simplicity and taking into account the analysis done in the previous section, in the 
following 6 = is taken in equations (p3|) and (0). 



A. Double-7 band heads 

There are three nuclei that have double-7 bands identified without ambiguity: ^^^Dy, 
^^^Er, and ^^^Er |[7|-[T0|. In this section the band heads of these nuclei are studied using the 
Hamiltonians (p3| - p^) to get an improved description of the anharmonicity phenomenon. 

The number of bosons for ^^^Dy and ^^^Er is = 16 while for ^^^Er it is = 15. In 
the different calculations shown in this section the parameters of the Hamiltonian have been 
chosen as to reproduce as well as possible not only the heads of single- and double-7 bands 
but also the structure of the bands. Because these calculations are schematic and do not try 
to be the best answer, the parameters of the Hamiltonian, for simplicity, will not be given 
fully and only in the case of final spectra will be shown. For comparison, the calculation 
with three-body terms (see section |ITl|and reference [|^) is also included. 
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In figures ^ H, and ^ the heads of single- and double-phonon bands are shown. In each 
figure six panels are included: a) experimental data, b) calculation with Hamiltonian ( p3D 
and X = ~a/7/2, c) calculation with Hamiltonian ( p3D and x = —0.55, d) calculation with 
Hamiltonian (p^) and x = — e) calculation with Hamiltonian (p^) and x = —0.55, 
and f) calculation with a three-body term (rf^ x x (i^)^^) ■ (d x d x d)'^^\ In each panel, 
from left to right are represented the ground state, the 7 band head, the (3 band head, the 
double-7 K = Q band head, and the double-7 K = A band head. Due to the controversy on 
the nature of the j3 band, several candidates for the latter have been included. For the same 
reason in the figures will be used the label . In ^^^Dy and ^^^Er no information on the 
double-7 K = band exists. The value x = —0.55 is chosen as an alternative to the SU{?)) 
value because it describes very well the E2 transition probabilities in this mass region. Also, 
the same value of x is taken in the Hamiltonian and in the transition operator, in line with 
the Consistent-Q Formalism (CQF) [^. In section |V C| this ansatz is used in the complete 
analysis of spectra and B{E2) transitions for ^^^Dy, ^^^Er, and ^^^Er. 

The most striking feature of figures is that in all calculations the position of the 
double-7 band heads and the degree of anharmonicity is well reproduced. Thus, the energies 
of the different band heads only are not sufficient to completely determine the Hamiltonian. 
Nevertheless, not all possible terms are able to create sufficient anharmonicity in the double- 
7 bands. For example, in the case of three-body terms, a phenomenological study of the 
most relevant type of term shows that only x x d'^)^'^^ ■ (dx dx d)'^^^ is able to produce 
the required anharmonicity (see section |T| and reference On the other hand, only a 

few four-body interactions have been explored here and it cannot be excluded that other 
four-body terms can produce the appropriate degree of anharmonicity. 

In order to decide which Hamiltonian is more appropriate, a description should be at- 
tempted not only of band heads but also of the structure of the bands and of E2 transition 
probabilities. 
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B. Moments of inertia 



The study of the moments of inertia of the lowest bands is a very sensitive way to test 
the different calculations shown in figures ^ and ^. Particular attention will be paid to 
the dynamic moment of inertia which can be obtained from the relation between angular 
momentum and 7-ray energy [^,^] and can be approximated by, 

T^2rf§., (27) 

where J is given dimensionless. In a plot of 7-ray energy versus J, 1 will be the slope. 
Equation (l27| ) can be used to study the structure of the rotational bands in comparison with 
experimental results. 

In figure |l^ the moments of inertia of the ground and 7-bands in the nuclei ^^^Dy, 
^^^Er, and ^^^Er are compared to those obtained with the Hamiltonian of section |T|. The 
parameters for ^^^Dy are k = 24.2 keV, x = —0.55, = —6.0 keV, and ^4 = 51.0 keV, for 
^^^Er are indicated in section ^ and for ^^^Er are k = 24.1 keV, x = —0.55, = —1.9 keV, 
and 64 = 31.6 keV. The predicted moments of inertia disagree completely with the almost 
pure rotational structure observed experimentally. So, although this simple description 
based on one single three-body term describes well the anharmonic position of the double- 
phonon band heads, it fails in the moments of inertia of ground and 7-bands. A more 
realistic description needs others three-body terms. A recent analysis shows, however, 
that Hamiltonians with two different three-body terms cannot get the correct moment of 
inertia. 

The results with the Hamiltonian p3|) are shown in figure |ll|. Two different calculations 
are shown which correspond to panels b and c, respectively, of figures 0-^. In this case a 
better agreement is obtained but still some discrepancies remain, especially when a realistic 
value for x is taken {x = —0.55). Finally, in figure |12| the results with the Hamiltonian ( p^ 
are compared with the data. Again, two different calculations are shown which correspond 
to panels d and e, respectively, of figures Here, good agreement is found both for 

X = -V7/2 and x = -0.55. 
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The different results with the Hamiltonians (^) and p^), can be understood quahtatively 
by analyzing the structure of {Q{x) ■ Q{x)){Q{x) ' Qix))- With just two-body terms the 
Hamiltonians ( |23D and (^) are equivalent. This is different when up to four-body terms are 
included because {Q{x) ' Q{x)){Q{x) 'Qix)) does not only contribute to C'2[5'f/(3)]^ and 
but also to (L2)2. This can be clarified with the equation 

(Qix) ■ QixmQix) ■ Qix)) = \c,[su{3)]l - p,[su{3)]lP + ^{Lr. (28) 

As a consequence, ■ Q{x)){Q{x) " Qix)) substantially modifies the rotational structure 

of a band even in the case of pure SU{3). This is how can be qualitatively understood that 
a description in terms of pseudo-Casimir operators is the most appropriate for dealing with 
anharmonic vibrations. In the next section a complete analysis is given of the nuclei under 
study in the framework of CQF using the Hamiltonian (0). 



C. Realistic calculations 

The complete calculated spectra of ^^^Dy, ^^^Er, and ^^^Er and the most relevant E2 
transition probabilities are presented in this section. They are compared with existing data. 
In each calculation the same value of x has been used both in the Hamiltonian (^) and in the 
electromagnetic operator (^21) . Finally, the effective charge in the transition operator (^2|), 



Cgg, was fixed for each nucleus to reproduce the B{E2] 2i Of) value. 

The experimental and calculated spectra for ^^^Dy, ^^^Er, and ^^^Er are shown in fig- 
ures |14|, and |15|, respectively. The parameters used in the calculations are listed in 
table |I|. The parameters that yield the best fit to the energy spectra are very similar in the 
three cases which is consistent with the analysis carried out in the preceding sections. The 
overall description of the energies is satisfactory. The calculated low-lying bands are in good 
agreement with their experimental counterparts (see also figure |1^ while the double-7 band- 
head energies are close to the experimental values. The calculated ratios (H) are: Rq = 3.31 
and RJ = 2.81 for ^^^Dy, R^ = 3.08 and Rj = 2.43 for ^^^Er, and R^ = 2-93 and Rj = 2.46 
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for ^^^Er, to be compared with the experimental ones: Rj = 2.84 for ^^^Dy, Rq = 2.82 and 
RJ = 2.45 for ^^<^Et, and Rj = 2.50 for '<^^Er. Only for the 7^ vibration the experimental 
and theoretical results are slightly different in the sense that this framework overestimates 
the anharmonic behavior for the Jx^o band. This can be corrected by increasing the value 
of Ixl in the Hamiltonian which, however, will introduce one more parameter because in the 
electromagnetic operator a different value of x must be used. 

For the calculation of E2 transition probabilities the same x values as in the Hamiltonian 
are adopted. The effective charges are e|g = (1.66)^ W.u., e|g = (1.83)^ W.u., and e|g = 
(1.67)2 W.U., for ^6^Dy, ^^^Er, and ^^^Er respectively. In tables |IT| 0, and the observed 
B{E2) values and ratios concerning 7-vibrational states are compared with the theoretical 
results. In general, a good overall agreement is obtained in the three cases under study. 



VI. CONCLUSIONS 

In this paper the problem of anharmonicity in the j3 and 7 vibrations of deformed nuclei 
was addressed in the context of the interacting boson model. The occurrence or not of 
anharmonicity was shown to be related to the order of the interactions between the bosons 
and the conclusions of the analysis can be summarized as follows. If the Hamiltonian includes 
up to two-body interactions, no sizeable anharmonicity can be obtained and the observed 
behavior cannot be obtained. The origin of this behavior is related to the existence of a 



first-order phase transition between rotational and vibrational nuclei which excludes 
any interplay between one- and two-body terms, necessary to obtain anharmonic spectra. 
For up to three-body interactions that preserve SU{3) symmetry it can be shown that 
the observed anharmonicity cannot be fully reproduced. Furthermore, extensive numerical 
calculations indicate that even a general IBM Hamiltonian that includes up to three-body 
interactions has difficulty in reproducing all observed aspects of ground, single- and double- 
7 bands. However, due to the large parameter space of three-body interactions which is 
difficult to search exhaustively, this cannot be considered as a firm conclusion and alternative 
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approaches {e.g. mean-field) should be tried to tackle the same problem. Finally, it was 
shown that a simple parametrization of the IBM Hamiltonian that includes up to four-body 
interactions can account for all observed properties of the three deformed nuclei with firmly 
established double-7 vibrations. In particular, the introduction of SU{3) pseudo-Casimir 
operators allows to describe the double-phonon states while keeping correct the properties 
of low-lying bands. 

It should be emphasized that, in spite of its fourth-order character, the Hamiltonian 
considered here is only slightly more complex than the usual IBM Hamiltonian (one more 
parameter) and is a straightforward extension of the Consistent-Q Formalism that was pre- 
viously successfully applied to many nuclei. Also clear from our study is the need for more 
experimental information about the double-phonon vibrations in deformed nuclei beyond the 
three cases known at present: our analysis indeed has shown that this information represents 
a challenging test of any theoretical description of deformed nuclei. 
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FIGURES 




2.0 I 
1 2 

e/K 



FIG. 1. The ratios R^^ (as defined in the text) as a function of 9i/k for different I. 
Hamiltonian (|9|) is used with x = —0.5; the boson number is = 15. The dashed lines give 
experimental values for the corresponding ratios in ^^^Er. 
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FIG. 2. Experimental (a) and calculated (b) spectrum for ^^^Er. The theoretical results are 
obtained with the Hamiltonian (^) with k = 23.8 keV, x = —0.55, k' = —1.9 keV, and 64 = 31.3 
keV. The boson number is = 15. 
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FIG. 5. The ratio R^^^^'^ (as defined in the text) for the Hamiltonian 

H = a C2[5C/(3)] + c C2[SU{?,)f with = 5, 10, and 15. 
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FIG. 6. Experimental (a) and theoretical (b) spectrum for ^^^Er. The theoretical results are 
obtained with the Hamiltonian H = 13.43 - 20.84 C2[SU{3)] + 9.296 IQ-^ C2[SU{3)]'^ (all 
coefficients in keV) with x = —V7/2. The boson number is N = 15. 
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FIG. 7. Band heads of ground, 7, /?, double-7 K = and double-7 K = 4 bands of ^^^Dy. 
Panels correspond to (a) experimental data, (b) calculation with Hamiltonian ( p3| ) and x = —V^ /2, 
(c) calculation with Hamiltonian ( p3| ) and x = —0.55, (d) calculation with Hamiltonian p^ ) and 
X = — \/7/2, (e) calculation with Hamiltonian ( p^ ) and x = —0.55, and (f) calculation with 
Hamiltonian (U). 
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FIG. 8. Band heads of ^^^Er. See caption of figure]^. 
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FIG. 9. Band heads of ^^^Er. See caption of figure]^. 
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100 200 300 400 500 600 700 

(J^J-2) (keV) 

FIG. 10. Total angular momentum J (dimensionless) of the initial state versus 7-ray energies 
for AJ = 2 transitions in ground and 7 bands, for ^^^Dy, ^^^Er, and ^^^Er. Full lines correspond 
to experimental data and long-dashed lines correspond to the calculated results obtained with the 
Hamiltonian (P). The corresponding parameters are given in the text. 
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FIG. 13. Experimental (a) and theoretical (b) spectrum for ^^^Dy. The Hamiltonian (|2^ ) is 
used with parameters k' = 12.18 keV, a = —92.90 keV, c = 0.05150 keV, and x = —0.55. The 
boson number is iV = 16. 
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FIG. 14. Experimental (a) and theoretical (b) spectrum for ^^^Er. The Hamiltonian (]2^ ) is 
used with parameters k' = 13.55 keV, a = —75.40 keV, c = 0.05286 keV, and x = —0.45. The 
boson number is iV = 15. 



34 



(a) 



3,* 
2,' 



5* 

4* 



2 

T K=4 



Exp. 
Er 



(b) 



6,* . 

4,* . 
2* 

0,*- 



5,*- 
3,*- 



P 



6/- 



V e;. 

2e* 

0/ 

2 4^* 

T K=0 ' .2 



T K=4 



Theo. 
Er 



FIG. 15. Experimental (a) and theoretical (b) spectrum for ^^^Er. The Hamiltonian (24) is 
used with parameters k' = 13.23 keV, a = —67.25 keV, c = 0.04080 keV, and x = —0.50. The 
boson number is iV = 16. 
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TABLES 

TABLE L Observed and calculated B(E2) values and ratios for i^egj. ^ schematic calculation 
using an SU{2,) Hamiltonian. The E2 operator is used with e^^ = (1.97)2 W.u. and x = -0.26. 

i?(E2) value or ratio 



Observed 



Calculated 



5(E2;2+ ^ 0^) (W.u. 

S(E2;4+ ^ 2+) (W.u. 

S(E2;2+ ^0^) (W.u. 

ij(E2;0+ ^2+) 

5(E2;2+ ^ 0^) 
ij(E2;4+ ^2+) 

5(E2;2| ^ 0^^) 



214 ± 10 
311 ± 10 " 
5.5 ±0.4 " 
3.8 ± 1.3 ^ (2.2+q;7 

1.3 ±0.4 ^ (0.9 ±0.3 



214 
302 
5.4 
2.5 

2.5 



From reference [30| 
From reference 0. 
^ From reference IHl. 



TABLE IL Parameters of the Hamiltonian (|2J) obtained in the best fit to spectra and B{E2) 
transitions in the nuclei ^^^Dy, ^^^Er, and ^^^Er. 



Nucleus 



k' (keV) 



a (keV) 



c (keV) 



X 



N 



166^^ 
i68Er 



12.18 
13.55 
13.23 



-82.90 
-75.40 
-67.25 



0.05150 
0.05286 
0.04080 



-0.55 
-0.45 
-0.50 



16 
15 
16 
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TABLE III. Observed and calculated B{E2) values and ratios for ^^^Dy. The E2 operator (M) 



is used with 


e^g = (1.66)2 W.u. and x = -0.55. 










i3(E2) value or ratio 








Observed 


Calculated 


r 1 / -1 — 1 / -\ ,---v —1— 

B{E2; 2J 


0^) (W.u.) 


209 ± 3 " 


209 


B(E2;4]^ ^ 


2+) (W.u.) 


272 ± 14 " 


298 


S(E2;2+ ^ 


0+) (W.u.) 


4.0 ±0.4 " 


3.9 


i?(E2;4+ - 
5(E2;2+ ^ 


0^) 


0.5 - 3.9 ^ 


3.3 



" From reference 
^ From reference 1711. 



TABLE IV. Observed and calculated B{E2) values and ratios for ^^'^Er. The E2 operator 
is used with e^g = (1.83)^ W.u. and x = -0.45. 

i?(E2) value or ratio 







Observed 


Calculated 


S(E2;2+ - 


> 0+) (W.u.) 


214 ± 10 " 


214 


B{E2;4f - 


. 2+) (W.u.) 


311 ± 10 


304 


5(E2;2+ - 


. 0+) (W.u.) 


5.5 ±0.4 " 


5.9 


i?(E2;0+ - 

5(E2;2| - 
i3(E2;4+ - 

i?(E2; 2^ - 


-2j) 
-0+) 


3.8 ± 1.3 ^ (2.2+q;7 '=) 
1.3 ±0.4 ^ (0.9 ±0.3 '=) 


1.8 
2.7 



" From reference |^| 
From reference |^ 
From reference Isll 
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TABLE V. Observed and calculated B(E2) values and ratios for ^^^Er. The E2 operator (22) 
is used with e^g = (1.67)^ W.u. and x = -0.50. 

B(E2) value or ratio 





Observed 


Calculated 


B{E2;2f of) (W.u.) 


207 ± 10 " 


207 


B{E2;4t 2f) (W.u.) 


318 ± 10 


294 


B(E2;2+ ^Ot) (W.u.) 


4.80 ±0.17 " 


4.6 


i?(E2;4+ ^2+) 


0.5 - 1.6 


2.7 


5(E2;2+ ^ of) 






i3(E2;5+ ^3+) 


0.7-3.5 


2.1 


5(E2;2j ^ 0^") 







^ From reference ||32|| . 
^ From reference ITO . 
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